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Abstract 

The red colors of many solid bodies in outer Solar System may be caused by 
tholins, which are refractory organic complexes, incorporated in their surface materials. 
Tholins synthesized in the laboratory are shown to match the colors of these bodies when 
their optical properties are used in rigorous scattering models. We review recent 
successes in modeling the spectra of icy outer Solar System bodies with tholins as the 
coloring agents. New work on the systematic laboratory synthesis and analysis of tholins 
made by cold plasma discharge in mixtures of gaseous CH 4 /N 2 shows that the 
composition of the tholin depends strongly on the pressure in the reaction chamber, and 
only weakly on the mixing fraction of CH 4 relative to N 2 . In tholins made at high 
pressure (e.g., 23 hPa) the abundance of aliphatic hydrocarbons is greater and the 
abundance of aromatic hydrocarbons is less than in tholins made at low pressure (e.g., 

0. 13 hPa). Tholins made at low deposition pressures show a greater abundance of N-H 
bonds. 

Keywords: 

1. The Colors of Small Bodies in the Outer Solar System 

Solid bodies in the outer Solar System (OSS) exhibit colors ranging from neutral 
to very red. Color is defined as the gradient in the spectral reflectance with increasing 
wavelength, after correction for the color of the illuminating source, the Sun. Objects of 
low albedo (0.01-0.15) commonly show ared slope to varying degrees and extending 
somewhat longer than 1 pm. Icy planetary satellites with higher albedo (0.2 - 1.0) also 
have varying degrees of color, but the positive spectral reflectance gradient of these 
bodies is usually limited to the wavelength region 0. 2-0.6 pm. The Moon is red because 
of a space weathering phenomenon in which particles in the solar wind sputter and reduce 
small amounts of iron in the iron-bearing surface minerals. The reduced nano-phase Fe is 
deposited on adjacent grains, altering the spectral reflectance and rendering it red in color 
(Pieters et al. 2000, Hapke 2001). The same space weathering phenomenon occurs on 
many asteroids in the Main Belt (Hapke 2001). Mars is red because of its surface coating 
of oxides of iron (e.g., Binder and Cruikshank 1966), and the red and orange colors of Io 
result from the various polymorphs of sulfur covering its volcanic regions. Both Mars 
and Io have much higher albedos than the objects we address in this paper. 
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2. Materials that May Impart Red Color ' 

While a modest degree of red color, such as that of Trojan asteroid 624 Hektor 
can be explained by the presence of the mineral Mg-rich pyroxene (Cruikshank et al. 
2001; Emery and Brown 2004), redder objects require other materials that are not 
minerals. Gradie and Veverka (1980) introduced the concept of organic solids to explain 
the red color of the D-type asteroids, suggesting the presence of “...very opaque, very 
red, polymer-type organic compounds, which are structurally similar to aromatic-type 
kerogen”. They compared the reflectance spectra of the D-type asteroids to mixtures of 
montmorillonite (clay), magnetite (iron oxide), carbon black, and a coal-tar residue 
(kerogen) that was insoluble in organic solvents, finding a satisfactory match for both the 
low albedo and the red color of this specific class of asteroids. Kerogen-like structures 
are found in carbonaceous meteorites (Kerridge et al. 1987), and the presence of similar 
structures in interstellar dust grains has also been deduced from infrared spectra of dusty 
diffuse interstellar clouds (Pendleton and Allamandola 2002). 

Tholins, which are the refractory residues from the irradiation of gases and ices 
containing hydrocarbons, have color properties that make them reasonable candidates for 
comparison to the spectra of Solar System bodies. Furthermore, they are produced in 
laboratory conditions that are reasonably analogous to the conditions of the exposure of 
atmospheric gases and surface ices in planetary settings in various natural environments, 
particularly in the outer Solar System. A number of tholins have been prepared in the 
context of the photochemical aerosols in Titan’s atmosphere (e.g., Khare et al. 1984, Coll 
et al. 1999, Ramiriz et al. 2002, Tran et al. 2003, Imanaka et al. 2004), and while they are 
optically (spectrally) quite similar to Titan, they have until recently proven difficult to 
analyze and characterize fully from a chemical point of view. The recent work, for 
example by Tran et al. (2003), and Imanaka et al. (2004), has greatly improved the 
analysis and characterization of a class of tholins produced by photolysis and cold plasma 
irradiation of gaseous mixtures, as discussed below. 

Optically, the tholins are characterized by strong absorption in the ultraviolet and 
visible spectral regions giving them strong yellow, orange, and red colors, high 
reflectance at longer wavelengths, and (in some cases) absorption bands characteristic of 
aliphatic and aromatic hydrocarbons with varying amounts of substituted nitrogen. We 
discuss this further in a later section. 

3. Models of Planetary Surfaces 

While simple spectral comparisons can be a guide to the presence or absence of 
tholins or other materials on a planetary surface or in an atmosphere, quantitative studies 
require the computation of the full spectral reflectance properties, which include not only 
the spectral shape, but the absolute reflectance (characterized by geometric albedo) across 
the full range of wavelengths for which observations exist. Such quantitative spectral 
matching is usually accomplished by the use of rigorous scattering models, of which the 
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Hapke (1981, 1993) and Shkuratov et ai. (1999) theories are often used. In addition to 
the original literature sources, the discussion of the practical application of Hapke theory 
to reflectance spectroscopy of solid surfaces by Verbiscer and Helfenstein (1998) is 
particularly useful. See also McEwen (1991). Poulet et al. (2002) have compared the 
Hapke and Shkuratov theories in the context of modeling solid surfaces, and Cruikshank 
et al. (2004) have described the computational procedures in more general terms. 


Planetary surfaces and atmospheres are observed remotely both by the sunlight 
that is scattered from them and by their thermal emission arising from their natural 
temperatures. The specific wavelength regions of these two regimes depend on the 
characteristics of the individual Solar System body and on its distance from the Sun 
(hence its temperature). Typically, however, the scattered sunlight regime extends from 
the ultraviolet to about 5 pm wavelength, and the thermal emission extends from about 5 
pm to 200 pm. In this paper we are particularly interested in the spectral region where 
scattered sunlight dominates the flux detected from a Solar System body. 

We note that the naturally occurring materials on a planetary surface can be 
segregated from one another or can be mixed in various combinations in a number of 
ways that affect the scattering properties and the abundances derived from the synthetic 
spectra calculated from scattering theories to match the observed spectra. These mixing 
scenarios have been discussed in various publications, including Poulet et al. (2002), 
Cruikshank et al. (1998), Roush and Cruikshank (2004), and Cruikshank et al. (2004). 


4. Tholins Used in Spectral Models 

Models using tholins as coloring agents successfully match the spectra of a 
number of both unusual and representative OSS bodies. Cruikshank et al. (1998b) 
matched the spectrum of the extremely red Centaur object 5145 Pholus with a model 
using the Khare et al. (1984) Titan tholin, H2O ice, CH3OH ice, and amorphous carbon; 
the Pholus spectrum has clear absorption bands in addition to the red color. Owen et al. 
(2001) calculated a simple model of the spectral reflectance of the low-albedo 
hemisphere of Saturn’s satellite Iapetus (0.4 - 3.8 pm) that successfully matches the red 
color and a strong absorption band (2.7-3. 6 pm) with the Triton tholin of McDonald et al. 
(1994) [the optical constants of Triton tholin, 0.2-10 pm are shown in Fig. 3 of 
Cruikshank et al. 2005]. Khare’s Triton tholin was made by corona discharge in a 
gaseous mixture of N2 + CH4 (1000:1) that approximately simulates the tenuous 
atmosphere of Neptune’s satellite Triton. Dotto et al. (2004) matched the spectra of 
several Centaurs and Kuiper Belt objects having a wide range of color from nearly 
neutral to extremely red, also using the Titan tholin of Khare et al. (1984) in various 
combinations with other components (H2O ice, amorphous carbon). 

In a study of the rings of Saturn, Cuzzi and Estrada (1998) found that the particles 
in the A and B rings contain a material imparting a red color (absorption toward the violet 
spectral region), while the particles in the C ring and Cassini’s division are lower in 
albedo and less red in color. They note that, “No silicates have the appropriate 
combination of steep spectral slope and high absorptivity to explain the rings’ visual 
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color while remaining compatible with microwave observations.” Titan tholin matches 
the colors and albedos of the particles when incorporated into the scattering models with 
H 2 O ice. To explain the lower albedo and more neutral color of the darker rings, Cuzzi 
and Estrada (1998) suggest that “material with properties like carbon black, as seen in at 
least some comets and interplanetary dust particles, is needed. . .”. Subsequent modeling 
of the rings by Poulet and Cuzzi (2002) and Poulet et al. (2003) also incorporated tholins 
and amorphous carbon to achieve fits to the observational data (0.3-4 pm). 

5. New Research on Tholins 

Tholins are refractory organic materials formed from reduced gas mixtures or ice 
mixtures with UV photolysis and/or charged particle irradiation (Sagan and Khare, 1979). 
In order to calculate scattering models of the spectra of OSS objects, the optical constants 
of tholins are required for comparison with the astronomical observations. One of the 
most successful measurements of the optical constants was done for Titan tholin by 
Khare et al. (1984). They produced Titan tholin from CH 4 /N 2 = 10/90 gas mixture under 
DC plasma irradiation at 0.2 hPa and measured the optical constants of Titan tholin from 
soft X-ray to microwave frequencies. The geometric albedo of Titan (its atmospheric 
haze) was explained fairly well using the derived optical constants of Titan tholin 
(McKay et al. 1989). Since then, some sets of optical constants (the complex refractive 
indices, n and k) were measured for additional hydrocarbon tholins produced by charged 
particle irradiation of pure CH4 gas and CH 4 /H 2 gas mixtures (Khare et al. 1987), for 
tholins produced by charged particle irradiation of H 2 O/C 2 H 6 = 6/1 ice mixtures at 77 K 
(Khare et al. 1993), for poly-HCN (Khare et al. 1994), and for kerogen and organic 
residue from the Murchison meteorite (Khare et al. 1990). Although those optical 
constants show wide variations, all of these organic materials show reddish colors with 
greater absorption at shorter wavelengths in the visible spectral region. [Note on pressure 
units: 1 hPa = 1 hectoPascal = 10 2 Pascal = 0.76 torr] 

In the early work the cause of the reddish color of tholins and the relationships 
between the optical proprieties and the chemical structure of tholins were not well 
understood. Although polycyclic aromatic hydrocarbons and some functional groups 
were suggested in Titan tholin (Khare et al., 1984; Sagan et al., 1993), other details of the 
chemical structure remain unclear. Many laboratory simulations have been conducted to 
investigate the optical and chemical properties of tholins (McDonald et al. 1994; McKay 
1996; Coll et al. 1999, Ramirez et al. 2002). McDonald et al. (1994) compared the 
chemical properties of Titan tholin and Triton tholin that were produced from CH4/N2 = 
10/90 and 0.1/99.9 gas mixtures respectively, and found distinct properties in these 
tholins with specific elemental compositions and chemical structures. McKay (1996) 
showed that optical properties vary with the tholins produced from different initial 
CH4/N2 ratios by comparing their transmittances in the ultraviolet- visible (UV/Vis) 
spectral ranges. Ramirez et al. (2002) measured the optical properties of tholin formed 
from a CH4/N2 gas mixture (2/98) at about 1 hPa, and found a ten times smaller value of 
k (the imaginary component of the complex refractive index) in the UV/Vis ranges than 
those measured by Khare et al. (1984). This difference reveals that the properties of 
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ihulin may depend on the experimental conditions, such as initial ratio of CH4/N2 and 
deposition pressure. 

/ 

Recently, Imanaka (2004) and Imanaka et al. (2004) systematically investigated a 
range of Titan tholins produced under various pressures in cold plasma. They found that 
the ratio of aromatic to aliphatic components monotonically decreased with increasing 
pressure and that the k-value in the UV/Vis ranges correspondingly decreased by an order 
of magnitude. They suggested that the optical properties of tholin in the UV/Vis ranges 
may strongly depend on the delocalized electrons from aromatic components embedded 
in a saturated carbon matrix. 

The absorption properties of tholin-like materials that show a reddish slope in the 
visible wavelength region merit further consideration. Tholins are a carbonaceous 
material containing C, H, N, and O, and while the exact structure is unknown, the 
fundamental building blocks are most likely composed of carbon, and probably nitrogen. 
Because of the different bond hybridization available (sp, sp 2 , sp 3 ), carbonaceous 
materials show a wide variety of optical properties. For example, diamond is composed 
of tetrahedral sp 3 -hybridized C-C bonding configurations (a bonds) and is transparent at 
visible wavelengths. Graphite is a trigonal sp 2 network that forms planar six-fold rings. 
The completely delocalized 7i electrons are responsible for the metal-like properties of 
graphite (its black color and absence of reddish color). In general, the absorption in the 
visible range depends on the D Cf and n- lT excitations of unsaturated bond structures. 

The more conjugation in a molecule, the more the absorption is displaced toward longer 
wavelength (e.g., Rao, 1975). Nitrogen atoms might help the delocalization because of 
the lone pair electrons. Tholin is a highly heterogeneous material containing both 
saturated and unsaturated bonds. Thus, the reddish slope of absorption in the visible 
range would depend on both abundance and size distributions of conjugated structures 
(either aromatics or alkenes). Since tholin is formed from simple molecules, the size 
distributions of conjugated structures would be controlled by its formation chemistry. 
Hydrogen is one of the main elements that break C=C bonds, localizing the □ electrons, 
and thus determines the size of delocalized □ electrons. Thus, the role of hydrogen in the 
tholin formation (or modification) process would appear to be important in controlling 
the reddish absorption in UV/VIS wavelengths (Imanaka 2004). 

In order to investigate Solar System bodies in the full range of wavelengths that 
are observed, it is important to obtain the optical constants of various tholins over a broad 
range of wavelengths in laboratory, including the near-infrared (NIR, -1-3 pm) and mid- 
infrared (MIR, -3-50 pm). In particular, the absorption properties of tholin in NIR range 
have not been well investigated because of the general weakness of any absorption in the 
continuum or in specific bands. Figure 1 shows the transmittance spectra of Titan tholins 
formed at five different pressures from CH4/N2 (=10/90) gas mixture by cold plasma 
irradiation, which were previously described in detail (Imanaka et al. 2004). The near 
infrared spectrum of Titan tholin formed at 0.26 hPa is also shown. The sample thickness 
for NIR region is about twenty five times greater than those of other samples in MIR 
region, so that the absorption peak intensities in NIR are much smaller than in MIR. 
Several overtones of C-H, N-H, and CON stretching bands and combination bands are 
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preliminarily identified in the NIR range. The reddish absorption slope extends from the 
visible range to about 1.8 pm. It would be important to compare simultaneously the 
reddish slope in the UV/VIS ranges and absorption bands in NIR/MIR ranges using the 
optical properties of various tholins with observations of OSS objects over a wide range 
of wavelength. 



Wavelength (pm) 

Figure 1. Transmission spectra of Titan tholins formed by cold plasma discharge in a 
CH 4 /N 2 = 10/90 gas mixture at five different pressures. The spectra have been offset 
from one another for clarity. The NIR spectra of Titan tholin formed at 0.26 hPa was 
obtained from a sample about twenty five times thicker than the others, so that the 
absorption peak intensities in the NIR are much smaller than in the MIR. The reddish 
absorption slope extends from the visible range to about 1.8 pm. 


In a systematic study of the methane-nitrogen tholins, Imanaka (2004) has also 
investigated the variation in spectral properties of tholins made at a uniform pressure and 
temperature but with variable amounts of CH4 in the mixture. Figure 2 shows 
transmittance spectra of tholins formed at 1.6 hPa, but from different ratios of CH4/N2 in 
the initial gas mixtures. The ratios of the areas of the aliphatic C-H bands (~3.4 pm) to 
N-H bands (~2. 8-3.4 pm) increase monotonically, but only slightly, with the increase of 
CH4 ratio in the initial CH4/N2 gas mixture. This variation in the ratio of C-H bands to N- 
H bands is smaller than that in tholins formed at different pressures (see Fig. 1). The 
absorption properties in the UV/Vis ranges do not change significantly, except for the 
pure CH 4 tholin. Thus, the optical properties of tholin formed from CH4/N2 gas mixtures 
are dominated by the deposition pressures rather than their initial mixing ratio used in this 
study. 
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Figure 2. Transmission spectra of tholins formed by cold plasma discharge at 1.6 hPa 
from CH 4 /N 2 gas mixtures with various initial ratios. Intensities have been offset. The 
absorption intensity of N-H bands around 3.2 pm increases monastically when the CH 4 
ratio in initial gas mixture decreases. The absorption intensity of C-H bands around 3.4 
pm has a maximum when the CH 4 ratio is ~10 - 25 % in the initial gas mixture. 

Additional properties of the Titan tholins produced in this investigation were • 
discussed in greater detail in Imanaka (2004) and Imanaka et al. (2004), and are briefly 
summarized here: 

1) Both the quantity and size of aromatic ring compounds in tholin increase with 
decreasing deposition pressure. 

2) The abundance of saturated C-H bonds is lower in tholins formed at low pressures, 
while N-H bonds become more abundant. 

3) More nitrogen is included in carbon networks at lower pressures in structures such as 
carbon diimides (-N=C=N-) and isocyanides (-C-NDC), while at higher deposition 
pressures nitrogen exists as a terminal group, such as nitrile (-CCN). 

4) Tholins formed at low pressures contain clusters of N-containing polycyclic aromatic 
compounds in a matrix of C and N branched chain networks which are connected tightly 
to one another with hydrogen bonding of N-H bonds. Tholin formed at high pressure (23 
hPa) consists of polymer-like branched chain structures terminated with -CH 3 , -NH 2 , and 
-CON with fewer aromatic compounds. 

5) The C/N ratio in tholins increases from a value between 1 .5 and 2 at lower pressures 
to 3 at high pressure (23 hPa). Thus, the N is incorporated into tholin more efficiently at 
low pressures. 

6 ) Reddish-brown tholin films formed at low pressures (0.13 and 0.26 hPa) show 
stronger UV/Vis absorptions than yellowish tholin films formed at higher pressures (1.6 
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and 23 hPa). The imaginary part of the complex refractive indices at UV-Vis wavelengths 
depends strongly on the delocalized □ electrons from aromatic and N-containing 
heteroaromatic compounds in tholins. 
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6. Conclusion 

Together with ices and minerals, organic solid material is emerging as a 
significant component of the suite of constituent materials of small Solar System bodies. 
In this paper we have summarized some of the astronomical evidence, particularly in 
terms of the colors that even small concentrations of complex organic solids impart to the 
surface of an object. We have focused on recent laboratory work on the synthesis and 
analysis of a tholin formed by energetic processing of a gas mixture chosen to simulate 
the composition of Titan’s atmosphere, at a number of different pressures chosen to 
simulate different pressure levels (hence altitudes) in that atmosphere. This systematic 
study, reported in detail by Imanaka (2004) and Imanaka et al. (2004) has begun to 
elucidate the details of several aspects of the formation pathways and the chemical 
structure of Titan tholin. Energetic processing of other gas mixtures and of condensed 
ices also produces complex refractory organic solids of probable relevance to outer Solar 
System bodies, and while preliminary experiments with these materials have produced 
promising results, much more work is needed. 
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Figure Captions: 

Figure 1. 

Transmission spectra of Titan tholins formed by cold plasma discharge in a CH4/N2 = 
10/90 gas mixture at five different pressures. The spectra have been offset from one 
another for clarity. The NIR spectra of Titan tholin formed at 0.26 hPa was obtained 
from a sample about twenty five times thicker than the others, so that the absorption peak 
intensities in the NIR are much smaller than in the MIR. The reddish absorption slope 
extends from the visible range to about 2 pm. 

Figure 2. 

Transmission spectra of tholins formed by cold plasma discharge at 1.6 hPa from CH4/N2 
gas mixtures with various initial ratios. Intensities have been offset. The absorption 
intensity of N-H bands around 3.2 pm increases monotonically when the CH4 ratio in 
initial gas mixture decreases. The absorption intensity of C-H bands around 3.4 pm has a 
maximum when the CH4 ratio is ~10 — 25 % in the initial gas mixture. 
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